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Open charm and beauty production in hadron reactions 

G.I. Lykasov^, V.V. Lyubushkin^ and V.A. Bednyakov^ 
Joint Institute for Nuclear Research, 141980, Dubna, Moscow region, Russia 

The production of charmed and beauty hadrons in proton-proton and proton-antiproton colhsions at high 
energies is analyzed within the modified quark-gluon string model (QGSM) including the internal motion of quarks 
in colliding hadrons. It is shown that using both the QGSM and NLO QCD one can describe the experimental 
data rather successfully in a wide region of transverse momenta. We also present some predictions for the 
future experiments on the beauty baryon production in pp collisions at LHC energies and on the charmed meson 
production in pp reactions at GSI energies. 
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1. Introduction 

Various approaches of perturbative QCD in- 
cluding the next-to-leading order calculations 
(NLO QCD) have been applied to analyze the 
heavy flavour particle production in hadron reac- 
tions at high energies. 

Such reactions are usually analyzed within the 
hard parton scattering model (HPSM) suggested 
in |ll2j . This model was applied to the charmed 
meson production both in proton-proton and 
meson-proton interactions at high energies, see 
for example [3]. The HPSM is significantly im- 
proved by applying the QCD parton approach 
|4|5j , see details in [6] and references therein. Un- 
fortunately the QCD approach including the NLO 
has some uncertainties related to the renormal- 
ization parameters especially at small transverse 
momenta -pt [S] . 

In this paper we study the charmed and beauty 
meson production within the QGSM [7] or the 
dual parton model (DPM) [8] in -p-p and -pp colli- 
sions at high energies based on the 1 jN expansion 
in QCD |9|10j . We show that this approach can 
be applied rather successfully at not very large 
values of pt- In addition, we investigate the open 
charm and beauty baryon production in pp col- 



lisions at LHC energies and very small pt within 
the QGSM to find new information on the sea 
charmed and beauty quark distributions in the 
proton. And at the end of the paper we use this 
approach to analyze the charmed meson produc- 
tion in the pp collision at not large energies be- 
cause the obtained results would be very interest- 
ing for the future experiments at the GSI (Darm- 
stadt) planned by the PANDA Collaboration. 



2. Charmed and beauty hadron produc- 
tion in pp collisions at high energies 

2.1. Heavy flavour meson production 

First, let us analyze the D and B meson pro- 
duction in the pp collisions within the QGSM in- 
cluding the transverse motion of quarks and di- 
quarks in colliding protons [TT . As is known, 
the cylinder type graphs for the pp collision pre- 
sented in Fig.l make the main contribution to 
this process [7]. The left diagram of Fig.l, 
the so-called one-cylinder graph, corresponds to 
the case where two colorless strings are formed 
between the quark/diquark {q/qq) and the di- 
quark/quark (qq/q) in colliding protons; then, af- 
ter their breakup, qq pairs are created and frag- 
mentated to a hadron, for example, D meson. 
The right diagram of Fig.l, the so-called mul- 
ticylinder graph, corresponds to creation of the 
same two colorless strings and many strings be- 
tween sea quarks/antiquarks q/q and sea anti- 
quarks/quarks q/q in the colliding protons. 
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Figure 1. The one-cylinder graph (left diagram) 
and the multi-cylinder graph (right diagram) for 
the inclusive pp hX process. 



The general form for the invariant inclusive 
hadron spectrum within the QGSM is |12lll] 



E 



da 



2E* da 



d'^p TTy/s dxdpl 



^o-„(s)(/)„(x,pt) , (1) 



where E, p are the energy and the three- 
momentum of the produced hadron h in the labo- 
ratory system (l.s.) of colliding protons, E* , s are 
the energy of h and the square of the initial en- 
ergy in the c.m.s oipp; x and pt are the Feynman 
variable and the transverse momentum of h; an is 
the cross section for production of the n-Pomeron 
chain (or 2n quark-antiquark strings) decaying 
into hadrons, calculated within the "eikonal ap- 
proximation" [13], the function 0„(a;,pt) has the 
following form [TT]: 



4>n{x,pt) = I dxi / dx2llJn{x,PuXi,X2) , (2) 



where 

-lpnix,Pt;Xi,X2) = 

F^^^ {x+ , Pt ; Xi)F^:^ (x_ ,pt;x2) /F^^^ (0, ) 
+ F^:^ {x+,pu x^)F(^^ (x_ , PuX2)/f!^^ (0, Pt) 
+ 2(n-l)F^(;)(x+,pt;xi) 

xF!^\x^,Pt;x2)lF^:\0,pt) , (3) 



and x± = {y/x'^ + x^ ± x)/2, 
xt = 2^{ml+pi)/s, 



F^^\x±,pt;x^^2) 

= J d^ktfi^\xi,2,kt)Gr^h(^^,kt;pt^ , 

(4) 

F^'^Ho^Pt) 

= f dx'd^ktf^^'>{x',kt)Gr^h{0,Pt) 

Jo 

= Gr^h{0,Pt). (5) 

Here r means the flavor of the valence 
(or sea) quark or diquark, fr'^\x',kt) is the 
quark distribution function depending on the 
longitudinal momentum fraction x' and the 
transverse momentum kt in the n-Pomeron 
chain; Gr^h{z,kt;pt) = zbr-,h{z,kt;pt), 
Dr^h{z,kt;pt) is the fragmentation function 
(FF) of a quark (antiquark) or diquark of flavour 
T into a hadron h (D meson in our case). We 
present the quark distributions and the FF in 
the factorized forms fr{x,kt) — fr{x)gr{kt), 
and, according to [15, Gr-,h{z,kt;pt) = 
GT^h{z)gT^h(kt)i where kf = pt - zk*. We 
take the quark distributions frix) and the FF 
Gr^h{z) obtained within the QGSM from [T2l 
116117] . whereas their kt distributions are chosen 
in the form suggested in [14115] (see the details in 
i) 

g^kt) = {B^j27r) exp(-B,A(fci)) , (6) 
-gr^hik) = {.Bl/2n) eM-BcA{kt)) , (7) 



where A{kt) — {y^kf + rnjj — mn)) and rn\n = 
pI + rnjj After the integration of eq. ^ over d'^kt 
we have, according to [T3] , 

F^"Hx±,pt;xi,2) 

^ fi"\xia)Gr^h{z)In{z,Pt) , (8) 

where z = x±/xi^2, In{z,pt) = B^/(27r(l + 
B^uid)) e-Kp{-B^{mDt - mo)), B^ = Bc/{1 + 
npz'^), p = Bc/Bq. The function Bz also can be 
presented in the equivalent form Bz = Bq/(p + 
nz"^), where p = Bq/Bc- The differential cross 
section da /dpi for D mesons produced in pp col- 
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lisions is written in the form I 12j 

da TT / \ f I , \dx , 

To describe the experimental data on the pt- 
spectra in the pt region where the NLO QCD cal- 
culation has a big uncertainty we chose p — 7 
both for the Tevatron and the LHC energies. 
When Be < I (GeV/c)-\ Eq.® can be ap- 
proximately presented in the form gr—>h{kt) — 

l(a} + kX) at k\ < 2tod, where a = ^2mBlBc 
[l9] . This form is similar to the form for the FF 
of heavy quarks obtained within the perturbative 
QCD. Note that the function /„(z,pt) in Eq.® 
was obtained in |15lll| on the assumption of the 
consequent sharing of the transverse momentum 
Pt in the proton (antiproton) between n-Pomeron 
chains. It allowed us to describe rather satisfac- 
torily the experimental data on the inclusive pt 
spectra of charmed and beauty mesons produced 
in pp collisions at moderate values of the trans- 
verse momentum < 10GeV/c[5]. As is shown 
in [T3], this version of the QGSM describes rather 
satisfactorily the experimental data on inclusive 
spectra of D mesons produced in pp collisions at 
a/s = 27.4 GeV. 

It allows us to make the predictions for inclu- 
sive Pt spectra of and mesons produced 
in the pp collision at LHC energies and compare 
our results with the NLO QCD calculation for the 
produced charmed quarks [H], see Figs. 2, 3. The 
solid lines in Figs. 2, 3 correspond to our calcula- 
tions within the QGSM, whereas the hatched re- 
gions show the calculations within the NLO QCD 
including uncertainties [H]. A big difference be- 
tween the QGSM and NLO QCD calculations at 
Pt > 10 GeV/c for D and B mesons can be due 
to the following. First, the NLO QCD calculation 
[18] does not include the hadronization of quarks 
to heavy mesons, whereas the QGSM calculation 
includes it. Second, we do not include the con- 
tribution of gluons and their hard scatterings off 
quarks and gluons which can be sizable at large 
values of pt . 

2.2. Heavy flavor baryon production 

Now let us analyze the charmed and beauty 
baryon production in the pp collision at LHC en- 



1 0^ 1 PP->D°X, s'°=14 TeV, a^=0.0 
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Figure 2. The inclusive spectrum for Z?" mesons 
produced in the pp collision at the LHC energy 
^ = 14 TeV obtained within the QGSM for 
charmed mesons [6] and the NLO QCD for c 
quarks [l8| . 



ergies and very small pt within the soft QCD, 
e.g., the QGSM. This study can be interesting 
for it may allow predictions for future LHC ex- 
periments like TOTEM and ATLAS and an op- 
portunity to find new information on the distri- 
bution of sea charmed (c) and beauty (b) quarks 
at very low Q^. According to the QGSM, the 
distribution of c(c) quarks in the nth Pomeron 
chain (Fig. 1 (right)) is, see for example [11 and 
references therein, 

X(l - a;)"''('')-2aB(0)-H(ap(0)-a^(0))+n-l ^ (^0) 

where (5c(c) is the weight of charmed pairs in the 
quark sea, C'c(e) normalization coefficient 

[12] . a^(0) is the intercept of the tp- Regge trajec- 
tory. Its value can be —2.18 assuming that this 
trajectory a^s,{t) is linear and the intercept and 
the slope (0) can be determined by drawing the 
trajectory through the ^f-meson mass = 3.1 
GeV and the x-meson mass m^. = 3.554 GeV [2U] . 
Assuming that the ip- Regge trajectory is nonlin- 
ear one can get a^{0) ~ 0, which follows from 
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Figure 3. The inclusive spectrum for mesons 
produced in the pp colhsion at the LHC energy 
7i = 14 TeV obtained within the QGSM for 
beauty mesons [5] and the NLO QCD for single b 
quarks [18]. 



Figure 4. The differential cross section da/dx for 
the inclusive process pp AcX at -y/s 62 GeV. 



perturbative QCD, as it was shown in |21j. The 
distribution of b{b) quarks in the nth Pomeron 
chain (Fig. 1 (right)) has the similar form 



X(l - a;)"''(")-2aB(0) + (ap(0)-aT(0))+n-l ^ (^l) 

where ap(0) — 1/2 is the well known intercept of 
the /9-trajectory; asiO) ~ —0.5 is the intercept 
of the baryon trajectory, ar{0)) is the intercept 
of the T- Rcgge trajectory, its value also has an 
uncertainty. Assuming its linearity one can get 
ar{0) = —8, —16, while for nonlinear (bb) Regge 
trajectory q;t(0) ~ 0, see details in il7j. Inserting 
these values to the form for /"(g) (a;) and f^^i^) 
we get the large sensitivity for the c and b sea 
quark distributions in the nth Pomeron chain. 
Note that the FFs also depend on the parameters 
of these Regge trajectories. Therefore, the knowl- 
edge of the intercepts and slopes of the heavy- 
meson Regge trajectories is very important for 
the theoretical analysis of open charm and beauty 
production in hadron processes. 

The information on the charmonium (cc) and 



botomonium (bb) Regge trajectories can be found 
from the experimental data on the charmed and 
beauty baryon production in pp collisions at high 
energies. For example. Fig. 4 illustrates the sen- 
sitivity of the inclusive spectrum da/dx of the 
produced charmed baryons Ac to different val- 
ues for a.^(0). The solid line corresponds to 
a^(0) = 0, whereas the dashed curve corresponds 
to a^{0) = —2.18. Unfortunately the experimen- 
tal data presented in Fig. 4 have big uncertainties; 
therefore, one can't extract the information on 
the a^{0) values from the existing experimental 
data. 

A high sensitivity of the inclusive spectrum 
da/dx of the produced beauty baryons Ab to dif- 
ferent values for ar{0) is presented in Fig. 5 (left). 
The pt-inclusive spectrum of Af, has much lower 
sensitivity to this quantity, according to the re- 
sults presented in Fig. 5 (right). Actually, our re- 
sults presented in Fig. 5 could be considered as 
some predictions for future experiments at LHC. 
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CU/ = 



daPP 

= (Ti(s)[(l - uj)<j>PP{x,pt)+uj<j>{x,pt)] 



n=2 



(12) 



Figure 5. The differential cross section da/dx 
(left) and da/dP^ (right) for the inclusive pro- 
cess pp AbX at y/s — 4 TeV. 
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3. Charmed and beauty hadron produc- 
tion in pp collisions 

3.1. Heavy flavour meson production at 
high energies 

The production of heavy mesons like D and 
B mesons in proton-antiproton collisions at high 
energies is usually analyzed within the different 
schemes of QCD. To study these processes within 
the QGSM we have to include at least one ad- 
ditional graph corresponding to the creation of 
three chains between quarks in the initial proton 
and antiquarks in the colliding antiproton, as is 
illustrated in Fig. 6(c). 

The diagrams in Fig.6(a,b) are similar to the 
one-cylinder and multicylinder diagrams for the 
pp collision in Fig.l with a following difference. 
In the pp collision two colorless strings between 
quark/diquark (q/qq) in the initial proton and 
antiquark/antidiquark {q/qq) are created. Many 
quark-antiquark [q — q) strings for pp collision 
(Fig. 6b) are the same as for thepp collision (Fig.l, 
right diagram). Therefore, the invariant inclusive 
spectrum of hadrons produced in the pp collision 
calculated within the QGSM has the following 
form: 



Figure 6. The one-cylinder graph (a) and the 
multi-cylinder graph (b), and the three-chains 
graph (c) for the pp — > hX inclusive process. 



where 1 — w is the probability of contribution 
of the cut one-cylinder (one-Pomeron exchange) 
and cut multicylinder (multiPomeron exchanges) 
graphs (the left and right diagrams in Fig. 6), 
whereas lo is the probability of the contribution 
of the three-chain diagram to the inclusive spec- 
trum. The value of w can be estimated as the 
ratio of the pp annihilation cross section CTp^" to 
the total pp cross section cr*^* . The cross section 
(7*^* is well known in the wide range of the initial 
energies to the Tevatron energy, whereas exper- 
imental data on (Tp|^" are available only for the 
antiproton initial energy up to 10 GeV, see [21] 
and references therein. However, some theoret- 
ical predictions, for example [25 26 , show that 
asymptotically Cp^" goes to about 2 — 4 mb. It 
corresponds to a; ~ cr°^"/crp^* < 0.1 at the Teva- 
tron energy Note that in addition to the graph of 
Fig. 6c there can be diagrams consisting of these 
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three chains and multicilynder chains between sea 
quarks and antiquarks. However, as our esti- 
mations show, their contribution to the inclusive 
spectrum is much smaher than the contribution 
from the three-chain graph (Fig. 6c). Therefore, 
we neglect it. 

for the function 
(t>n{x,pt) entering 



The form 
is similar to 
by replacing 



into (3) 

(")/ 



(a;_,Pt;x2), i^i, (0,pt) 
to F^^\x_,pt;x2), F^q\Q,pt) respectively. 



and replacing F^q\x_,pt;x2), F^q'{0,pt) to 



in), 



F^"\x - , pt] X2) and ^g"''(0,pt) respectively. The 
additional term (l){x,pt) in (fT2)) has the following 
form 



4>{x,pt) = 3Fg^{x+,pt)Fg^{x.,pt)/Fg{0,pt) , (13) 



= F 



(n=l) 



{x±,pt) and 



where Fg^(^g^){x±,pt) 

Fg{Q,p,)=Ft='\0,Pt). 

The inclusive pt spectra of and mesons 
produced in the pp collision at the Tevatron en- 
ergy ^/s — 1.96 TeV are presented in Figs. (7,8), 
see [6]. The hatched regions in Figs. (7.8) show 
the calculations within the NLO QCD including 
uncertainties [29]. Note that the our calculation 
showed that the contribution of the three-chain 
graph (Fg.6c) is very small at the Tevatron en- 
ergy. It is due to small values of the pp annihila- 
tion cross section at very high energies |25|26) . 

3.2. Charmed meson production at inter- 
mediate energies 

Let us now use the QGSM to analyze the inclu- 
sive charmed meson production in pp collisions at 
not large energies, less than 15-20 GeV because it 
would be interesting for future experiments at the 
GSI (Darmstadt) planned by the PANDA Col- 
laboration. At these energies the contribution of 
the three-chain graph (Fig. 6c) can be sizable be- 
cause the cross section of the pp annihilation is 
not small '24j. Note that at energies close to the 
threshold of the D-meson production the binary 
process pp DD is dominant, according to [30] . 
Therefore, the total cross section of D mesons 
produced in pp collisions at energies starting from 
the threshold is the sum of cross sections for the 
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Figure 7. The inclusive P(-spectrum for 
mesons produced in the pp collision at the Teva- 
tron energy ^/s — 1.96 TeV obtained within the 
QGSM ^ (the solid hne) and within the NLO 
QCD [29 (the hatched regions). The experimen- 
tal data are taken from [27128] 



binary process Up^^i)]^ and the inclusive reaction 



pp^DX ■ 



pp- 

The cross section (Jpp- 



■DD 



was calcu- 



lated in [3TJ| within the simple Regge pole model. 
In Fig. 9 the cross sections (Jpp^DD (dashed curve) 
and (Jpp^DX (dash-dotted line and dotted line) 
are presented as functions of the incident mo- 
mentum of the antiproton in the l.s. of pp. The 
dash-dotted curve corresponds to our calculation 
of the inclusive process pp D{D)X with ne- 
glect of the graph of Fig. 6c (w = 0), and the dot- 
ted line corresponds to the similar calculation but 
with the Fig. 6c graph included; the probability 
for the three-chain graph (Fig. 6c) uj = cr?p"/cr^°* 
was taken from the fit of the experimental data 



Open charm and beauty production in hadron reactions 



7 



p p -> B X, s^"= 1.96 TeV, a = 0.0 



QGSM,tt) = 0.1 

■ Tevatron CDFII 
Uncertainty NLO 




5 10 15 20 25 
p,, GeV/c 



zi. 



10 



10 



10 



pp > DD 

■ pp > DX.m^O 

pp > DX,a = a(PiJ 

■ total cross section, Oi^ (£i(pj^ 



10 



50 



Figure 9. The cross section for D{D) meson pro- 
duction m. p — p collision as a function of the an- 
tiproton momentum piab in the l.s. 



Figure 8. The inclusive pt-spectrum for 
mesons produced in the pp collision at the Teva- 
tron energy ^/s — 1.96 TeV obtained within the 
QGSM (the sohd line) and within the NLO 
QCD [29] (the hatched regions). The experimen- 
tal data are taken from [27l28j 



[M| . The solid line is the total yield of D mesons 
produced in pp collision. 

4. Conclusion 

We have shown that the modified QGSM in- 
cluding the intrinsic longitudinal and transverse 
motion of quarks (antiquarks) and diquarks in 
colliding protons allowed us to describe rather 
satisfactorily the existing experimental data on 
inclusive spectra of heavy flavour mesons pro- 
duced in pp collisions and to make some predic- 
tions for similar spectra at LHC energies. 

To verify whether these predictions can be re- 



liable or not we apply the QGSM to the analy- 
sis of charmed and beauty meson production in 
proton- antiproton collisions at Tevatron energies 
including graphs like those in Fig. 6c correspond- 
ing to annihilation of quarks and antiquarks in 
colliding p and p, and production of heavy flavour 
mesons. We got a satisfactory QGSM description 
{pt < 10 GeV/c) of the experimental data on 
pt spectra of Z?*^ and mesons produced in the 
pp collisions which were obtained by the CDFII 
Collaboration at the Tevatron [55]. 

To describe these spectra and make some pre- 
dictions for the future LHC experiments in a wide 
region of transverse momenta one can combine 
the "soft QCD" (the QGSM) at smaU values of 
Pt with the NLO QCD at large pt- 

At the Tevatron energies the contribution of 
the three-chain graph (Fig. 6c) to the inclusive 
spectra of heavy mesons is too small, as was 
shown in [5]; therefore, it can be neglected. How- 
ever, at the antiproton energies about a few GeV 
the three-chain graph contribution is sizable and 
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can amount to 30-40 percent, as is shown in Fig. 9. 
It has also been shown that at the incident mo- 
menta of antiprotons piab above 8 — 9 GeV/c 
the inclusive production of D mesons in pp col- 
lisions should be included in addition to the bi- 
nary pp — > DD process to get the total yield of 
D mesons at piab < 14—15 GeV/c. These results 
would be interesting for future experiments at the 
GSI (Darmstadt) with the antiproton beam. 

We also made some predictions for future LHC 
forward experiments on the beauty baryon pro- 
duction in pp collisions which can give us new 
information on the beauty quark distribution in 
the proton and very interesting information on 
the Regge trajectories of (bb) mesons. 
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